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Thermodynamic Relations in n-Variable 
Systems in Jacobian Form: 


Part II, Polyphase Polycomponent 
Chemical Systems. 


F. H. CrawForp 
Williams College 


1. INTRODUCTION: THE FUNDAMENTAL JACOBIAN. 


In Part I* a general analytical theory of thermodynamic rela- 
tions was developed and applied to a variety of physical systems. 
In the present paper the same general technique is applied to a 
variety of chemical systems. 

Since it is possible to state the general methods and results in 
a form somewhat more suitable for the present case, a brief state- 
ment of the general procedure is necessary. 

The starting point is the expression for dU, the differential of 
the internal energy. This may be written in the form 


dU == + eee F de, (1) 


where the notation is the conventional one in thermodynamics. 
Thus T and S are, respectively, the absolute temperature and the 
entropy, while F;, ... F,, are the generalized forces applied to the 
system and 6;,... 0, are the conjugate geometric variables. Since 
in the present instance the term F,d6; will be —pdV, the hydro- 
static work term, it is convenient to rewrite Eqn (1) in this form 


dU = TdS — pdV + F,d6, +... + Fi dé, 
or l (2) 
dU = TdS — pdV + SF, d0,1 =n —2 
4=1 


Here / = n — 2 is simply the number of differential terms (if any) 


1F,.H., Crawford, This Journal, vol. 78, pp. 165-184, 1950. See also Phys. Rev. 
72, 521A, 1947, and Am. Jour. of Physics 17, pp. 1-5, 1949. 
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above the number required for the simplest hydrostatic case. 
Finally m is the total number of degrees of freedom or the 
number of independent variables required to describe the status 
of the system completely. Consequently Eqn (2) must always 
be assumed reduced to the minimum number of terms com- 
patible with the nature of the system. 

As in Part I the variables occurring on the right side of Eqn (2) 
are the primary variables for the problem, with the various thermo- 
dynamic potentials such as U, the enthalpy H, Helmholtz’s free 
energy A and Gibbs’ potential, G, as secondary variables. In 
addition a number of subsidiary variables will be introduced as 
convenience dictates. 

While in principle there is a wide choice of independent vari- 
able sets, for the most part it is convenient to make two restrictions: 

(1) Each independent variable will be one of the 27 primaries. 

(2) No conjugate pairs will be included in the set (i.., if T 

is chosen then S must not be, etc.). 
In the sequel a suitable set of independent variables (21, 2, ... Xn) 
will be called an x-set and unless otherwise stated will be a non- 
conjugate primary set. Once the «-set is decided upon, the re- 
maining primaries constitute the dependent variables (Yi, Yo, 
or simply the Y-set. 

In the selection of the non-conjugate x-set one variable alone 
must be taken from each differential term in Eqn (2). If the first 
variable is taken from each term (i.e., T, p, Fi, ...F,) then these are 
said to be selected in regular order. Had (T, p, Fi, 02, F3 ...) 
been selected, then T, p, Fi, Fs, etc., are selected in regular order 
while 62 is selected in reverse order. For reasons that will become 
apparent shortly, it is more convenient to have regular order 
determined by the order of the x-variable selected from some 
special term called the key-term. Thus if F2d@2 were taken as 
the key-term and 62 taken as x2, this selection would by definition 
fix regular order and the remaining T, p, Fi, F3, etc., would be 
taken in reverse order. 

It was shown in Part I that a certain Jacobian, the so-called 
fundamental Jacobian, J,, could, once the x- and Y-sets were 
selected, be made the core of the subsequent treatment. This 
Jacobian is by definition simply the Jacobian of the Y’s with 
respect to the x’s, where the subscripts are in the order of the 
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terms of Eqn (2) from which the variables come, counting from 
left to right. Thus 

_ 9 (V1, Yo, ... Yn) 
0 (x1, oe Ma) 


Jo (3) 


and contains 7? partial derivatives. Depending on the signs oc- 
curring in Eqn (2) and the precise way the x-set is chosen, we 
have some pairs of elements, symmetric about the main diagonal 
of the determinant, which are equal. Since these equalities are 
important in reducing the number of quantities in J, which must 
be measured (directly or indirectly) experimentally, it is very 
important to be able to render the entire Jacobian symmetric. It 
becomes therefore of importance to modify J, in such a way as 
to insure the symmetry for each choice of x-set. This new Jacobian 
(which may be identical with J, and will always be so if all the 
terms of Eqn (1) are positive and the x-set is all selected in the 
same order) will be called the symmetric or modified Jacobian, 
J’., for the problem. Like J,, it consists of a square array of n? 
partial derivatives. Since, however, it is symmetric we can by 
inspection write down equations expressing the equality of deriva- 
tives symmetrically placed about the main diagonal (i.e., conjugate 
pairs in the square array). These equalities are m(m-+1)/2 in 
number and represent the independent Maxwellian relations or 
the equations for the “ reciprocal effects” characteristic of the 
system. 

‘To insure J’, being symmetric and, as a matter of practical con- 
venience, to avoid as many minus signs as possible, we proceed 
as follows. Since we can always choose all of the x-set, or all 
of the Y-set, as positive, let us agree to take the x-set positive and 
place any necessary minus signs in front of the conjugate Y’s. 
It is convenient to give the steps for modifying J, in a simpler 
way than those given in Part I and at the same time to apply them 
to a specific example. Thus given dU = TdS — pdV +- F2dé2: to 
setup J’,. First: Select the desired x-set, say (T, V, 62). Second: 


0 (S, p, Fe) 
d(T, V, 62) 


Set up J., Jo = . Third: Select a key-term. This 


is often TdS but may be any other positive term. It is selected so 
as to introduce as few minus signs in J, as possible (using the sign 
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rules below). Let us then take TdS as the key term. This at once 
establishes regular and reverse order. Thus T is the first of the 
pair (T,S) from the key-term, whereas both V and 62 are the 
second members of their parent terms. Hence both V and 62 are 
selected in reverse order. Fourth: Modify J, by the following 
sign rules for the Y-variables. Each Y is given either (a) the 
sign of its parent term (when the conjugate x-variable was selected 
in regular order) or (b) the reversed sign of its parent term (when 
the corresponding x was selected in reversed order). 


(S, p, —Fe2) 
d(T, V, 62) 


only symmetric but has only one column of minus terms. Had we, 
on the other hand, taken F2d62 as the key term we should have had 


Thus we have J’, = , a Jacobian which is not 


oS, = ®) , which is symmetric to be sure, but contains twice 
(T, V, 02) 


as many minus signs. 


2. RESTRICTED SYSTEMS. 


In Part I the argument was carried out on the assumption that, 
in principle, any 7 of the 27 primaries could be taken for an x-set 
(the so-called unrestricted system). When, for various reasons, 
this is not so we have a restricted system and the specific restric- 
tions must be kept clearly in mind. Thus if F; in Eqn (2) should 
be a function of T alone it could obviously never be in the same 
x-set with T. In general any variable which depends functionally 
on less than 7 other variables, is a variable of restricted functional 
dependence (a restricted variable for short). Restricted variables 
are usually chosen as Y-variables and thus automatically excluded 
from the x-set. 

Although no a priori reasons exist for expecting a preponder- 
ance of one type of system over the other in nature, a casual 
survey of known systems renders it probable that restricted sys- 
tems are much more numerous than unrestricted. Thus, systems 
involving surface tension, reversible cells, temperature radiation, 
and many dielectric and magnetic systems, as well as general 
chemical systems, all fall in the restricted class. It would appear 
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then that the general unrestricted system is the exception rather 
than the rule. 

Returning to Eqn (z) let us select an x-set which contains no 
restricted variables, say, as an example (T, V, 6;, 62,...6:). Then 
the same x-set would be unchanged in form, as would J’,, even if 
the restriction should later be removed. 

We have with F;d6, as the key-term: 


—dS —dp oF, OF; 

a 

—dS —dp dF, OF; 

CO 

_ 0 (-S, —p, Fi,... Fi) —dp OF see (4) 

(T, V, 61, 6;) 06, 06; 06; 06; 
—dp OF; OF, 

| 06, 00, 00, 


(Note that had we taken TdS as the key term we should have had 
only two columns of plus terms.) Because of the symmetry of 
J’. Maxwell’s relations can be written down simply by inspection, 
beginning with 


p should be a restricted variable depending, say, on 


(T, 611), then (0p/00;) TV0, 
would of necessity vanish, as would its conjugate, and we should 
have two zeros occurring in J’,. Due to the symmetry any off- 
diagonal zero automatically introduces another symmetrically 
placed, whereas zeros on the main diagonal may occur singly 
and have no effect on the number of Maxwell’s relations. In the 
unlikely event that p depended on T alone we should have (7—1) 
zeros in the second column. Since the first of these would be on the 
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main diagonal the remaining (7—2) would be off diagonal and thus 
(n—2z) Maxwellian relations would become trivial. Consequently 


n(n 


instead of having — n(n — 1)/2 = derivatives to 


determine from experimental data, we should have only 


2 


and instead of 2(m — 1)/2 non-trivial and independent Maxwellian 
relations only m(m — 3)/2 + 2. 


3. SINGLE Pure SUBSTANCE OF VARIABLE Mass. 


Let us now examine these systems in which the masses of the 
various constituents are to be regarded as variable. Before pro- 
ceeding to the general case of chemical systems it is first of ad- 
vantage to examine the special case of a single pure substance 
under hydrostatic pressure. Here Eqn (2) becomes 


dU = TdS — pdV + pdm (5) 


where by definition » = (dU/dm), ,, and is the electrochemical 
potential or the Gibbs potential of a unit mass of the material. 
Forming the partial derivative of U with respect to m (T, p con- 
stant) gives us 

(6) 


=p=—u-—Ts — pu 


T'p Tp Tp 
where u, s and v are the corresponding quantities per unit mass 
or the specific internal energy, entropy and volume, respectively. 
Since p is defined then in such a way as to depend on specific and 
intensity variables only, it is independent of volumes and masses 
and therefore constitutes a restricted variable and we have here 
a simple example of a restricted system. 

Although we are at liberty to choose any legitimate set of 
independent variables as the x-set, let us choose, say, (T, p, 7). 
The Y-set is then (S, V, ») and with pdm as the key term we 
write out 


< 
€ 
= 
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—0S OV 
oT oT oT 
J’ d (—S, Vip) Op 
"A(T, p,m) | dp ap ap 
—0S OV Op 
om om am 


mov 
—mos 
=| dp op (7) 
—S Y O 


Since p = p(T, p), (02/8) rp must vanish. The other equivalents 
in the second determinant come from the relations du = Tds — pdv 
and du = —sdT + vdp and the definitions s = S/m, v = V/m and 
u = U/m. As a consequence we have only 5 derivatives which 
must be measured: 0s/dT, dv/dT, dv/dp, s and v, while two of 


the Maxwellian relations reduce to identities. 


This somewhat 


surprising reduction of two relations to identities must be ascribed 
to the rather special definition of » rather than to the general 


behavior of restricted systems. 


In contrast to the above simple situation, if we choose as the 
x-set (IT, V, m), J’, becomes more involved. 


Thus 
aS ap 
«OT 


J’ 0(S, —p) aS op 


° | aVe ave 
aS ap 
Om om 


(8) 
E 
dv’ 


mos op 


where in the second determinant the variable set has been changed 
from (T, V, m) to (T, v, m). This is readily done by the 
Jacobian method as an example will show. Thus take (0S/dm) 1, . 


(S,V) (ms, mv) 
_ O(m,v)r 9(M,V)r _ 
0 (mM, v)r 0 (mM, v)r 
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We notice that, although no element of this new modified Jacobian 
vanishes we still require only five partial derivatives (since s 
and v are really (0S/dm) 7) and (dV /dm) 7p respectively) and have 
only one significant Maxwellian relation. This illustrates the fact 
that although care in the selection of the x-set may reduce the 
complexity of results, it cannot reduce the actual number of 
physical parameters which must be determined experimentally. 


4. MULTICOMPONENT CHEMICAL SYSTEMS. 


In chenwcal systems with more than one chemical component 
in equilibrium we have restricted systems having interesting varia- 
tions of the properties discussed above. Now while this discus- 
sion does not produce any new results or results not deducible 
by conventional means, it does offer them in a compact and easily 
deduced form which is notably lacking in the usual presentation.’ 
The Jacobian methods are further very flexible and enable changes 
of variable to be made so as to minimize intermediate steps and 
the consequent likelihood of error. 

In a system of c chemical components (where c is by definition 
the number of chemical species whose masses can be varied in- 
dependently), the c chemical potentials p3, po, ... ue are always 
restricted variables. Depending on the number of phases present, 
other variables such as p or even T become restricted and the 
selection of the particular x-set which will be of greatest service 
in a particular case becomes of importance. We must therefore 
examine the situation in detail even at the expense of repeating 
much that is familiar to the chemical reader. 

In the first place we must distinguish very clearly between 7, 
the number of intensity variables which are independent, and 
n, the total number of independent variables of all sorts. Here 
m, <n and, for chemical systems under no other external force 
variable than pressure, is given by Gibbs’ phase rule 


(where c is the number of chemical components and ¢ the number 
of distinct physical phases present). 


2See for example R. W. Goranson, “ Thermodynamic Relations in Multicom- 
ponent Systems,” Carnegie Institute of Washington, D. C., 1930. 
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Useful as this rule is it does not fix 7 nor does it give us help 
in selecting x-sets in which less than the maximum number of 
intensity variables is utilized. In the case of the one component 
single phase system discussed in Sec. 3 we clearly have n = 3 and 
n,;— 2. The first x-set, (7,p,72), included the maximum number 
of intensity variables (2) and one extensive variable (m7). In 
the second x-set, (T,V,72), one intensity and two extensive vari- 
ables appear. But what happens when we go to a two-phase 
system? The value of 7; obviously decreases to unity, but does 
nm remain the same? And what happens when we have two or 
more components? 


5. VARIABLES AND CONSTRAINTS IN CHEMICAL SYSTEMS. 


From a somewhat superficial standpoint one might argue that 
a discussion of this sort is tautological. For in Eqn (2) of Sec. 1 
if we set ] = c, the number of independently variable masses, and 
consider T and p as two extra variables, then clearly 7 = / + 2, 
since it was explicitly stated that the number of terms in Eqns (1) 
and (2) was supposedly reduced to 7, the number of independ- 
ent variables required! But if this is so, what happens when p 
becomes a restricted variable? Does the number of terms in 
Eqns (1) and (2) depend on the number of phases present or not? 
Clearly we must consider the question of restrictions on the 
masses of the chemical species in our system as well as the method 
of determining the value of c in a given case. 

We may regard the question of mass restrictions on a specific 
collection of chemicals under two headings: (a) variable total 
mass, and (b) constant total mass. The second case can further 
be broken up into (bi) constant mass systems with semipermeable 
walls and so arranged that on withdrawing any one species we 
compensate by adding equivalent masses of one or more other 
species, and (bz) isolated systems where the mass is constant be- 
cause the walls are unpermeable to all chemical species. If we 
examine the definition of the number of components, c, naively, 
we might be tempted to say that if, in the case of a given set of 
non-reacting substances, the value were c for case (a) it would 
by definition be c-1 for case (b,) and zero for case (be), since 
here we have no direct way of varying amy mass at all. 
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It appears, however, that Gibbs regarded c as being fixed once 
for all as the maximum number of variable masses in the system. 
It is thus the maximum of potentially variable masses whatever 
later mass restrictions are imposed. This is clear from Gibbs’ 
original discussion of the relation n; = c—¢-+ 2.5 He repeatedly 
states that c — ¢ + 2 gives the greatest conceivable number of 
variations in phase of which the system is capable. ‘That is to say that 
1; = dmar — ¢, and the system can vary from 1 through 2, 3, 
etc., UP tO ¢maz phases, where this is given (when nm; = o), by 
mae — € + 2, and this value can be reached for a system without 
any reference to mass restrictions at all. 

The determination of the value of c in many complicated sys- 
tems is a difficult problem, though in the majority of cases met in 
practice rather readily settled. In the case of a collection of § dif- 
ferent stable and mutually inactive substances, c = §. In the case 
of systems with one or more chemical reactions present we may 
extend our name chemical substance to include any chemical spe- 
cies (molecule or fragment) present in the system. If we then 
have a number, 7,, of independent equations connecting these 
species (and these may include reactions of all sorts, ionizations, 
statements of electrical neutrality, etc.) we have c = 8 — n,. Just 
which particular molecule or fragment to choose as a component 
in a given case is a matter of taste and experience. Once a given 
species has been selected as a component one caution is necessary 
regarding computing its total mass and concentration in a given 
phase. Thus if in a gaseous system we take Hz as a component and 
the element existed at the same time as H and H:O vapor, we 
should have to add the total equivalent in both these forms to the 
free molecular Hz to get the total mass (and later, mass concentra- 
tion) of this component. 

With this in mind we may conveniently treat the system as 
though it were made up of c = § — n, non-reacting stable chemi- 
cal substances. 


6. DEGREES OF FREEDOM OF A CHEMICAL SYSTEM. 
Although Gibbs’ phase rule is invaluable in all such discussions 
we must find 7 itself as well. The process will incidentally indi- 


%J. W. Gibbs, Scientific Papers, Vol. I, page 76, et seq. (Longmans, Green and 
Company, 1906). 
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cate the subsidiary equations which must be used in all thermo- 
dynamic calculations as well as facilitate the selection of possible 
x-sets when these are not to be entirely primaries. 

The general procedure is simply to list the number of variables, 
Nvar, Which we need for the description of systems of ¢ com- 
ponents and ¢.phases. We then write down all the simultaneous 
and independent equations involving these variables. ‘These will 
include arbitrary definitions and external constraints, as well as 
experimental laws. Since we need only the total number, 74, of 
these equations and, in the case of the laws of nature, their gen- 
eral functional form, this listing can be done in a very general 
way. Then the number of degrees of freedom, , of the system is 
simply 2 = Nvar — Neg: 

The variables may conveniently be divided into the extensive 
variables (masses and volumes) and the intensive variables (con- 
centrations, pressures and all quantities which are functions of 
these alone). A collection of variables suitable for a c-component, 
o-phase system is given in Table I below. Here subscripts refer 
to components and primes to phases. Thus 773”, c2’, and ws” refer 
to the mass, mass concentration and chemical potential, respect- 
ively, of the second component in the third phase. When no prime 
occurs the property is that of the subscript component, where no 
subscript is present to that of the whole phase, and when neither 
is present to a property of the whole system. In the third column 
of the table is given the total number of variables of the group in 
question. 
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Table I 
Group Variables Number 
I , , 
111, 1119, Me, V 
(Extensive) 
I 
ms... me m* V* 
Mle, Me m, V 
II 
(Intensive) 
(c + 2) 
, 
Ill 
(Intensive) 
1s +++ 
, , 
eee He 
Total 3¢c+ 29+ 4 


Adding the quantities in the right column gives us 
Nvar = 3¢C + 26 + € + 4 


as the total number of variables with which we are concerned. 
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Table II 
Group Equations Number 
Mass and 
volume sums (j = 1, 2, $) 


j 


= m, 2 
j j 


In the same way in Table II are listed the independent equations 
connecting these variables, with the number of each type in the 
third column. The sum total of these gives us 2¢¢ = 3¢¢ + 24 + 2. 
Consequently* we have for 7 

N = — Neg = C+ 2 (12) 
and the value of 7 is independent of the number of phases present. 


We are therefore justified in rewriting Eqn (2) as a generalized 
form of Eqn (5) with / = c terms of the form pdm, i.e., 


dU = TdS — pdV + >> wy am, (13) 
4=1 
with the knowledge that 7 and hence the number of terms is in- 
dependent of whether we have a single phase or the maximum 
number, @dmae = 2 = c + 2. 
On combining the phase rule with Eqn (12) we have 
m=c+2=—-%"+ ¢ 
or (intensive) (extensive) (14) 
*If any component is missing from a phase then that concentration and mass 


drop out of nvar while one equation each drops out of the last two groups of 
Table II. Consequently this difference remains constant. 


Equations of state Vi = (T, p, mi,...m), ¢ (10) 
Definitions of c’s cj = mi / m’, (i,j) 
Nature of p’s pi = w (T, ci, cee » (4, 9) $c 
Gibbs’ Criterion , | (11) 
Total 3¢¢ +- 26 + 2 
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which simply means that if we employ the maximum number of 
intensive variables we shall then need in addition a number of 
extensive ones equal to the number of phases present. 

Since the above discussion assumed the total mass variable we 
now add the case (b) situation with 7z = constant and have 


case (b) ma md @- (15) 
as the general relations we sought. It is thus possible to take an 
x-set suitable for a case (a) system and by dropping out a suitable 
extensive variable (say 7 or V) have a suitable x-set for case (b). 
As soon as we pass to systems with two phases or more, primary 
x-sets cease to be as useful as those involving certain of the sub- 
sidiary variables from ‘Table II. 


7. RELATIONS IN SINGLE PHASE SYSTEMS. 


Although Eqn (13) is of the same form whatever the number 
of phases, due to the different choices of variables (and the con- 
sequent differences in subsidiary equations such as those in Table II 
and others) it is convenient to consider systems of single, double, 
and more, phases separately. 

To begin with let us write down J’, from Eqn (13) with 
primary x-set (T, p, 71, M72, ...™,). We have, with as 
the key term: . 


On Ope Ope 
of of or oF oT 
Ope 
op op op op op 
—oS oV Ope Ope 
a V, Itc) Om, Om Om} Om, O71, 
 0(T, p, m1, Ma, Mc) OV Ope Ope 
0M, OMe OM, O12 OM» 
OM, OM, OMe Om, 


(16) 
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a square array of (c + 2)* terms with something less than 
(c + 2) (¢ + 3)/2 independent. The (¢ + 2) (¢ + 1)/2 Max- 
wellian relations are read off at once by inspection. In cases where 
we want a few of the Maxwellian relations without the bother 
of writing out the square array we may do so from the differential 
symbol of J’, itself. Thus the main diagonal elements are given 
by forming the derivative of a Y-variable with the « immediately 
below it (such as —dS/dT, dV /dp, etc.). All the others are con- 
nected in pairs, where any desired pair is obtained by simply 
selecting the desired Y’s together with the x’s below them and 
equating the “cross” differentials. Thus let (—S, y1) be the 


desired Y’s. We write the scheme: and the wanted 


relation is simply —(0S/dn1;) Tp... = (0p1/0T) pm, 


Again with V and pe we have at once the array aa 7 5 and so 
6 


The procedure is satisfactory as long as the Y’s are taken with 
their proper signs, whatever their relative position in the differ- 
ential symbol for J’, may be. 

Another set of variables may be used in which the maximum per- 
mitted number of intensity variables, c +- 1, and a single extensive 
variable are employed. For this set we take (T,p,C1,Co, . . 
and eliminate the component masses 771,772, . . . #- from Eqn (13). 
Using the defining equations for the c’s from Table II together with 


m = um and des = 1 we have, after regrouping terms, 


(17) 


| | Pc) C1 + (p2 Pic) C2 Co_1 +- po 


| 
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or 
(18) 
dU = Tds — pdV + Bide, + Bodcot ... + Bos dec1+B.dm 


where the B’s are simply abbreviations for the coefficients of 
Eqn (17). We notice in particular that 


B; B,(7,p,C1,C2 Co_1,1) 
B. — B.(T,p,c,ce2 Coa) 


which therefore makes B, a restricted variable since it is inde- 
pendent of 7. 

We may now write out the new modified Jacobian in terms of 
this new x-set, where since B, is a restricted variable a zero occurs 
in the lower right-hand corner. We notice further that since 
0B,/0m = B;/m and by symmetry 0B;/dm = the lat- 
ter derivative is simply B,/m, etc., with 6B,/dm = o as the last 
element of the column. As a result the new Jacobian has a single 
zero on the main diagonal and (c-1) derivatives above this, which 
do not need to be measured as derivatives. 
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Although the Jacobians of Eqns (16) and (19) can now be 
used to calculate any desired thermodynamic property of the 
general chemical system, little information is gained which can- 
not more readily be got from specific examples. This is all the 
more true since in practice we seldom have more than two com- 
ponents and in such cases can readily use the appropriate part of 
the general results above. 


8. Two CoMPoNENT SINGLE PHASE SYSTEMS. 


If for purposes of illustrating the directness and flexibility of 
Jacobian methods we confine ourselves to systems of two com- 
ponents Eqn (13) reduces to four terms, since 7 = 4, and we have 

dU = TdS — pdV + pdm, + p2dMe (20) 
and with the primary x-set (T,p,71,m2) the modified Jacobian 
is a special case of that in Eqn (16). Hence 


OV Ope 
—dS dV a 
| OP Op Op) 
Om, OM, OM, 
—0S OV du, Ope 


Here the a’s are simply convenient abbreviations for the deriva- 
tives which, due to the symmetry, are reduced from 16 to 10 
independent ones. Each of the derivatives in turn represents 
some physical property of the system. Thus the top row gives 
temperature changes, the left column heat properties (essentially), 


a3 


433 


etc. In particular a) is written with a minus sign since then a) = 


(0S/dT) and T(0S/dT) as the heat 


capacity of the system at constant (,771,7712), usually a positive 
quantity. Of course many of the terms, in fact the four in the 


(21) 


| 
| 
4 
| 
| 
a 
¥ | 
| 
=. 
if 
= 
ees 
a 
13 
423 
& 
a 
. 
§ 
§ 
aN 


210 F. H. CRAWFORD 


upper left corner, can be written in other ways. Thus with 


S=ms and V=m1v, = 
s+ » (OV Tn, 
= v + (dv/dm) Tym’ etc. 


To determine any property involving other x-sets in terms of 
the a’s we proceed with the usual Jacobian transformation, to 
change from the new x-set to the old one, (7T,p,7711,772). If we 
take C Von’ the heat capacity at constant (V,72,772), we must 

1/442 


first notice that 


aT 


= Psp) 


Vm 


and we must essentially transform an entropy derivative in terms 
of the x-set (T,V,721,m72) to the set (7,p,711,7712). We have 


8(S,V) 
a(S, V 11,2) 
( ) (S, V 711,12) 0 ( 
oT 0 ( V ( V ) 
( M2) op Tm me 


where, as usual, variables written outside a bracket or a Jacobian 
are constant for all differentiations involved. We now need only 
read off the equivalents in the upper corner of the determinant of 
the a’s in Eqn (21) to have 


2 
T T (a0411 + a7) 
411 


(22) 


a result independent of the number of components. 
Again, to obtain the heat of addition per formula mass of, say, 
component 1 at constant (T,V,7722), we proceed as follows. If Mi 
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is the molecular “ weight ” of this component then dm, = M,dv,, 
if dy is the change in the number of moles of 1 corresponding to a 
change of mass dm. 


Thus TVm _ 1 TVm _ (08 
dm, M, dvy 0m, TVme 
0(S,T,V 
Om, TV me 0(m,T, V M712) V m2) 
0 (11, 
0(S,V) 
aV 
(=) 
p 


where the minus sign outside the numerator arises from inter- 
changing the order of (7,p) to give correspondence with that 
employed in Eqn (21). This gives finally 


dQ TVmo (41412 — 42431) 
411 (23) 


The chemical potentials »; and p2 are functions of intensity 
variables alone and hence, although depending directly on T and 
p, vary with 7, and mz only as these alter the concentrations. 
As a result the present x-set does not reveal the true character of 
the derivatives of the type (du;/d77 , (Opu2/dm , etc. 

Ivativ type Tom, (Ope/ Tym, etc 


For this purpose the second x-set, (T,p,c,7), used in Eqn (19) is 
preferable and we accordingly express these derivatives in terms 
of it.° The first of these gives us 


5In this section only the symbol c stands for c,, the concentration of the first 
component, and is not to be confused with c, the number of components which is 
here two, and of course not a variable. 
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— d(c,m,T,p) €,711) rp 

0(C,M) rp 


where the order of the new x-set has been written (c,77,T,p) 

simply to facilitate “ cancellation” without any troublesome sign 
changes. With #2; = mc and 22 = m(1—c) this result reduces to 


(1—C) (Op1/dc) , and likewise 


71 


Tpmye Tpm 


(9p2/81) = (1-—¢) (9p2/9C) 


om, 
(Op2/0m2) Tpm = (Opu2/dC) T pm 
—C 
(Op1/ 0712) Tpm, (Opu;/dC) Tpm 


Then since the second and fourth of these are equal by symmetry 
we need measure only one derivative, say (0u1/0C) 7pm (and not 
three). For this reason it is possible to deduce all the properties 
of this system from 8 rather than 10 fundamental derivatives 
say 41, 42, 43, 412, 413, and (0p,/0C) rpm). 

Another result is of interest. Since constant (7721, 7772) and con- 
stant (Cc, 7) are essentially equivalent as far as all derivatives in the 
two top rows of the Jacobian of Eqn (21) are concerned, all these 
derivatives are the same for the sets (7p,721,7772) and (T,p,c,m). 

The above discussion suggests that we can with profit turn to 
the form of Eqn (18) corresponding to two components, viz.: 


dU = TdS — pdV + Bide + Bodm. (24) 
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Here the modified Jacobian is 

(25) 
—0S 
oT oT by be dbz 


—0S OV OB, by Diy hig 


_ a(-S,V,B,,B2) | Op dp ap 


J, = 
d(T,p,cym) —aS 6B, @Bs| | Bt 
0c 0c 0c 0c mM 
bz; big 


072 Om Om 


Again examination shows that, as before, we need measure only 
8 actual derivatives for all the quantities in the square arrays to 
be determined. 

Since derivatives involving secondary quantities can be 
treated as readily as other derivatives, take the case of (0U/0S) yre 
and let us express it in terms of the b’s. As before 


0(U,V,T,c) a(U,V) 
dU 0(U,V,T,c) O(pm,T,c) re 

0(p,m,T,c) re 


Here primary derivatives come directly from Eqn (25) while 
partials of U must be got from Eqn (24) and then replaced by 
their equivalent b’s. The result is finally 


Drom big — hep — 


—bybi3 + 


0 


(Th, + + (Ths + phis — Be)bis 
— 


The present set of variables is well suited to the special case 
of fixed mass, for in this case we simply drop the last term 
of Eqn (24) and suppress the bottom row and right column of 


J’o, Eqn (25). 
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9. CoMPONENTs AND Two PuHasEs. 


The chief differences when multiphase systems are encountered 
arise from the increase in the number of restrictions with the con- 
sequent alteration in the possible x-sets. 

When we pass from one to two phases the set (T,p,7721,72, . . . Mc) 
no longer is possible, since p = p(T,¢1,C2,...Cc1). Since 
nm, = c — + 2 = c we must then omit either T or p (usually the 
latter) and introduce an extensive variable (often V). In the case 
of two components and (T,V 771,72) as the x-set we select 
as the key term and have 


—dp Op. Ope 
of of or of 


—dp Om, Ope 


—0S —dp Om. Ope 
0m, Om, Om, 


—dp Om Ope 
02 OMe 


Now unfortunately although we can read Maxwell’s relations 
directly out of Eqn (26), this set of variables is not as convenient 
experimentally as other sets made up chiefly of non-primary 
variables. A very convenient one for this purpose is (T’,ci,7',m") 
where c; is the concentration of the first component in phase one 
and mm’ and m” are the respective phase masses. For this change we 
need the following equations from Table II, as well as relations 
for S and V in terms of specific quantities. 


Mm S=wms + m's’ 


Then we may readily determine any derivative in Eqn (26) 
directly in terms of our new variables. 


oS 
Thus to express (7) , for example, in terms of 


(T,ci,m',m",) we have 
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0 (S, V ) 
( os ) 0 (S, V 11,M2) a ( T jsv 
0( 7 


where J8” is simply an abbreviation for the fourth order Jacobian 
in the numerator and Jif represents the third order Jacobian re- 
sulting from suppressing the first row and first column of JS’. 
We now write out JSY 


oS 0M, OMe 
or of or oT 


os 07714 
dci acl dg de dy de 


jr" = = (27) 
oS OV om 1 om 2 Cc 
om om om onl 


d, dz ds dg 


OS OV dm, ol 


om’ ont om’ ont" 


where dj, dz, etc., can be expressed in terms of masses and specific 
variables if desired. 


os’ os” 
(Thus MR + aT’ de = —ds, de = —dz, etc.) 


In the same way the other elements of the first column of the 
array in Eqn (26) can now be found. 
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The elements of the second, third, and fourth rows of the 
Jacobian of Eqn (26) need Jacobians which we shall designate as 
Je¥, Je” and J” respectively. Writing these out we have 


jev == V 


J 


and 


Juv 


(pr, V 721,712) 


O(p2, V 711,112) 
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They are seen to differ among themselves and from J*” only in the 
first column. Hence J;; is the same, whichever of the four Egns., 


(27) to (30), we use. Thus 


0( p,7711,M2) 
TmmMe 0( V 711,712) Ju 
0(ci,m' jm") 
0(p,V me) 
7722) 
Ope Ope 


Other useful sets such as (T, ci, x m) may be used where 
x = m"/m, the mass ratio in the second phase. Since this offers 
no special divergencies from the usual technique, no example 
need be worked out. 


10. Two CoMPONENTS AND THREE PHASES. 


When another phase is added to our system, 2; = ¢ — 3 +2 =I. 
Thus only one intensity variable is allowed and this we shall take 
as T. This gives the same variable set as used in Eqn (26), Le., 
(T,V 71,712), but here we have more rigorously restricted de- 
pendent variables. Thus p = p(T),u1 = pa(T),p2 = pe(T) and 
our modified Jacobian becomes much simpler. 


oy 
% | 
, 
4 
s 


218 F. H. CRAWFORD 


(31) 
—oS —dp Ope dp dus dps 
ot oT oT of oT oT oT oT 
y' d(—S,—p,p1 p2) _ 
° A(T, | —dS —dp du. Ape 
0m, Ot, OM, OM, 0m, 
—0S —dp Op. Ope 


Here we have 3 independent Maxwell relations and thus only 
4 independent derivatives in terms of which all other properties 
can be computed. 

As before, a more practical experimental x-set is given by 
(T, m',m’,m'”). Any of the derivatives in the first column of 
Eqn (31) can be got as ratios readily expressed in terms of a 
Jacobian of (S,V,7, m2) with respect to this new set. Thus 


a(S, ) where J; is the minor 

( oS ) _ _ Joo of J oo obtained by sup- 

oT Ji. pressing the first row 
Vmyme 

0( and column as before. 


Writing Joo out for reference gives 
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Here ci, cz and c3’ are included for symmetry even though 
cj + Ce = 1, etc., for the other two pairs. Further, the e’s are 
abbreviations for sums over the three phases, i.e., 


det and 


i 


Here of course the s’s, v’s and c’s are all functions of T only. 
Using the Jacobian Joo we can obtain the other S derivatives: 


0(S,711,M2) 
(m',m'm'") 
re) (S, V mz) 
Om, O(m, V me) OMe TVm Jus 


The same Jacobian Joo is also of use for changing from, say, the 
x-set (T,V,m',m") to (T,m',m" Thus 


with similar results when we start with (T,V,m',m'") or 

Other choices such as (T, x’, x”, m), where x’ and x” are the 
mass fractions in the first and second phases respectively, may be 
similarly treated. 
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We note finally that on suppressing the bottom row and right 
column of the array in Eqn (31) we have the appropriate results 
for the two phase single component system. 

In the case of two component systems with four phases we 
have invariancy, and T, 1, pe and all specific entropies and vol- 
umes are constant. Thus the only variables in Eqn (13) are 
S, V, mm, and 72, and these can be altered only through changing 
the masses of one or more phases. Consequently any desired 
quantities can be got more or less by inspection and although 
interesting transformative Jacobians can be set up, their use is 
rather artificial. 

It is hoped that the above treatment will serve to show the 
great directness and power of the Jacobian method of treating 
the thermodynamic properties of chemical systems. The results, 
further, involve a minimum of intermediate steps and auxiliary 
equations and, left in determinant form, are as compact and easy 
to handle as the complexity of the situation permits. 


F. H. CrawForp 
February 1, 1954 
Malaga, Spain 
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